Non-equilibrium phonon dynamics in trapped ion systems 
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We propose a concrete experiment to probe the non-equilibrium local dynamics of the one- 
dimensional Bose-Hubbard model using a trapped ion system consisting of a linear chain of few 
Ba + ions prepared in a state of transverse motional mode which corresponds to a fixed number of 
phonons per ion. These phonons are well-known to be described by an effective Bose-Hubbard model. 
We propose a protocol which leads to a sudden local sign reversal of the on-site interaction strength 
of this Hubbard model at one of the sites and demonstrate that the subsequent non-equilibrium 
dynamics of the model can be experimentally probed by measuring the time-dependent phonon 
number in a specific motional state of the Ba + ions. We back our experimental proposal with exact 
numerical calculation of the dynamics of a Bose-Hubbard model subsequent to a local quench. 

PACS numbers: 37.10.Jk, 37.10.Rs, 37.10.Ty, 63.20.kg, 64.70.Tg 



The possibility of the use of systems of ultracold 
atoms and trapped ions to emulate theoretical models of 
strongly correlated condensed matter systems has been 
intensely investigated in recent years [lHj]. One of the 
first examples in this regard constituted the realization 
of the Bose-Hubbard model in optical lattice systems and 
characterization of the superfluid-insulator transition of 
the model via measurement of momentum distribution of 
the bosons [l[ . More recently, there has been several the- 
oretical proposals of emulating various condensed matter 
models in trapped ion systems • The main advantage 
of the trapped ion emulators over their ultracold atom 
counterparts is that they have relatively large separation 
between individual ions compared to that between lattice 
sites of cold atoms in an optical lattice This allows 
optical measurement of the properties of the vibrational 
states of individual ions and thus enables one to address 
the local properties of the underlying strongly correlated 
model emulated by these systems. Also, these systems 
provide unique opportunity to study the Bose-Hubbard 
model in presence of an attractive on-site interaction. 

In recent years, it has also been realized that ultra- 
cold atomic systems, which are probably the best exam- 
ples of experimental realization of closed quantum sys- 
tems, provide an unique opportunity for studying the 
non-equilibrium dynamics of their constituent atoms near 
a quantum critical point ■ This has led to a plethora 



of theoretical studies of such dynamics However, 
in all of these experiments, it has not been possible to 
achieve protocols which takes the system out of equilib- 
rium by a local change of system parameters. This is 
mainly due to several experimental difficulties associated 
with addressing individual atoms in an optical lattice 
setup. Consequently, in most of the theoretical studies 
of non-equilibrium dynamics undertaken on these sys- 
tems so far [12l4l4| , the effect of a global time-dependent 
ramp of a system parameter has been studied in details; 
no attempts have been made to study the response of 



the system to, for example, a local ramp. Trapped ion 
emulators provide us an ideal test bed for studying such 
dynamics. 

In this work, we demonstrate, via proposal of a con- 
crete experiment, that a system of one-dimensional (ID) 
trapped ions provide a requisite experimental setup to 
study the non-equilibrium dynamics of a finite-sized ID 
Bose-Hubbard model subsequent to a local quench. The 
key differences between our proposal and the experimen- 
tal or theoretical studies carried out on cold-atom sys- 
tem earlier are two fold. First, our proposed setup pro- 
vides an implementation of a local quench which changes 
the on-site interaction of the emulated model at a given 
ion site (or equivalently at a lattice site for phonons). 
Second, it allows optical measurement of local physical 
quantities such as the time variation of projection of the 
phonon (boson) number into a motional state at a given 
site and hence provide us with direct information regard- 
ing local phonon (boson) number variation during non- 
equilibrium dynamics subsequent to the quench. We pro- 
vide a detailed description of the proposed experimental 
setup and chart out the parameter space in which such 
an experiment can be carried out without violating the 
phonon number conserving approximation. We also back 
our experimental proposal by an exact numerical study 
of the non-equilibrium dynamics of the finite-size Bose- 
Hubbard model emulated by the trapped ions. We com- 
pute the time evolution of the local boson (phonon) num- 
ber density n, subsequent to a sudden local sign change 
of the interaction (t/, — > — Ui), show that the variation 
of rii can be observed experimentally by measuring the 
time dependence of the motional ground state phonon 
occupation number at that site (n^o), and demonstrate 
that the amplitude of the quantum oscillations of rijo 
is maximal when J /Ui lies in the critical (crossover for 
finite-size system) region. We note that dynamical prop- 
erties of the Bose-Hubbard model emulated by an ionic 
system subsequent to a local quench has not been stud- 
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ied either experimentally or theoretically so far; our pro- 
posed experiment and the supporting theoretical analysis 
therefore constitute a significant extension of our under- 
standing of local non-equilibrium dynamics of correlated 
boson systems emulated by trapped ions. 

The experimental setup which we propose in the 
present work constitutes a linear chain of few Barium 
ions as shown schematically in middle and right panels 
of Fig. [TJ The relevant energy level diagram of a singly 
charged barium ion is shown in the left panel of Fig. 
[TJ We note here that most of our observations hold re- 
gardless of the nature of the used ions as long as they can 
undergo side-band cooling. For preparing the experimen- 
tal setup, we trap these ions in a linear trap operated at 
15 MHz radio-frequency and with the trap stability pa- 
rameter q ~ 0.42 used for the radial confinement. This 
will generate a pseudo-potential in the radial plane cor- 
responding to ion oscillation frequency lu x ~ 2.25 MHz. 
The confinement in the axial direction is only by DC 
voltages applied to the endcap electrodes. This can be 
made shallow so that the axial frequency is ~ 180 kHz 
and the mean distance between the ions can be made of 
the order to 20 /im 15]. These parameters will lead to a 
tunneling strength of J ~ 0.55 kHz. This in turn leads 
to fi x — 2J/uj x , the ratio between the coulomb interac- 
tion and the trapping potential, to be ~ 5 x 10~ 4 . In 
a linear trap like the one presented here, the inter ionic 
distance in the Coulomb crystal varies along the chain; 
however, for the present purpose this has been neglected 
as X <C 1. 

Next, we consider the on-site interaction for the 
phonons. For this, we note that for Doppler cooling of 
the ions, we use a diode laser at 493 nm. In addition, the 
side band cooling can be performed by driving the red- 
sideband near the S-D transition frequency at 2051 nm. 
In order to generate the on-site interaction, we use an off- 
resonant standing wave laser at 300 nm which results in a 
shift of the radial frequency uj x . This leads to on site in- 
teraction term U = 2(— l) s Fr] x , where F is the strength 
of the dipole interaction provided by the standing wave, 
and S = 0(1) corresponds to ions being trapped at the 
maxima (minima) of the trapping potential. For F ~ lo x 
and rf x = 5.26 x 10" 5 , we find a typical U = 235 Hz 
which leads to J/U ~ 2 . Note that since J, U -C uj x , the 
number of phonons are conserved. We stress that it is 
possible to tune J/U to as low as 0.09 by suitably choos- 
ing F ~ 25u> x with 300 nm standing wave and the radial 
frequency shifts to uj x — » 0.52uj x (1.31wa:) for repulsive 
(attractive) interactions; the system still remains in the 
phonon number conserving regime where Frj x << lo x . 
The standing wave pattern will be made by on-to-one 
focussing on each ion which eventually allows individual 
addressing. 

To implement a local quench on the prepared chain of 
Doppler cooled ions, we prepare a series of steps as shown 
in Fig. [5] The ion system, after being side band cooled to 
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FIG. 1: (Color online) (a) Relevant atomic levels of Ba + ions. 
The Zeeman sublevels are not shown for clarity, (b) Schematic 
representation of the trapped ion chain, (c) The harmonic 
oscillator levels in the radial trap. 

the motional ground state, is first prepared for simulation 
[step(a)]. A Raman type PI pulse with one blue sideband 
and one carrier photon applied to all ions then produces 
a state with one phonon each at all the sites [step (b)]. It 
is well-known, that the low-energy Hamiltonian of these 
phonons can be represented by the Bose-Hubbard model 
given by @ 

H = J^(b\bj +h.c.) + U^n i (h i -l)+u; x ^2h i) 

(ij) 1 * 

(1) 

where bi denotes the annihilation operator of the bosons 
(phonons) at site i and = b\bi is the local density op- 
erator. Note that Eq. [JJis identical to the Bose-Hubbard 
Hamiltonian in optical lattices with {7/2 — > U @. To 
emulate the effect of the quench of U at a given site, 
the applied single-site interaction produced by the stand- 
ing wave is suddenly quenched by changing the phase of 
the wave by means of a peizo actuator with frequency 
wo ~ 50 — lOOKHz. The frequency u>o is chosen such that 
it keeps the system in the phonon number conserving 
regime (ujo <C u> x ) and that it is fast compared to typical 
energy scales of the underlying Bose-Hubbard Hamilto- 
nian (w 3> J,U). This leads to the implementation of a 
local quench protocol. We stress that the ion-trap system 
provides a large enough window where quench dynam- 
ics of its underlying Bose-Hubbard Hamiltonian can be 
studied accurately within the number conserving approx- 
imation. In the rest of this paper, we restrict ourselves 
to this regime where uj x ^> uj ^> J, U. 

After the implementation of the local-quench protocol, 
we aim to measure the subsequent phonon distribution 
as a function of time. To this end, we note that the time 
dependence of the phonon distribution is reflected in the 
time-dependent local population of the phonons rijo(i) at 
a specified site (i) in a specified motional state, namely 
the ground state denoted by the index in this case. 
Thus we aim to measure the time-dependent population 
of the state \Si ,0) after the quench has been performed. 
In order to measure this population, a red-sideband PI 
pulse between Si and fls states are applied so as to 
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FIG. 2: (Color online) A sequence of laser pulses that are needed to be driven in order to implement site specific quenching, 
(a) Sideband cooling laser applied continuously in order to prepare the ID lattice in the ground state of the system, (b) A 
Raman PI pulse is applied in order to prepare the lattice in an initial near Mott-Insulator state. This pulse is applied to all 
the sites, (c) This is the final step after the quench pulse has been applied as described in the text. In this step shelving of all 
other states apart from the motional ground state is performed by a red sideband (RSB) PI pulse. 



shelve all other motional state population except that of 
|jSi,0). Note that the state |5i,0) is decoupled from 
this transition. Therefore a S — P fluorescence count will 
give the population in the motional ground state [step 
(c)]. The photon count rates can be measured within a 
time window of 1 second which is much larger than the 
life time of the excited state Pi/ 2 (7.8 ns). The total 
emitted photon count can be expressed as 

Ri(t) = (R)n i0 (t) =n l0 {t)fojrQ e Q /2, (2) 

where (R) is the mean fluorescence photon count rate; 
/ = 0.73 is the branching ratio between P\/2 and Si/ 2 
state, w = 0/4tt = 1/2 [1 - y/l - NA 2 } = 0.04 (where 
NA = 0.4 is the numerical aperture of a fluorescence 
collection Halo lens) is the fractional solid angle of de- 
tection, Q e is the quantum efficiency of a photomultiplier 
tube (PMT) detector which can be estimated to be about 
50% at wavelength of 493 nm and Q — 0.1 is the overall 
loss factor in collec ting the fluorescence photons in the 
experimental setup |16l |. Such a fluorescence count mea- 
surement as function of time therefore yields the requisite 
information on the time-dependent non-equilibrium pop- 
ulation distribution of the phonons subsequent to a local 
quench. In what follows, we shall provide a theoretical 
estimate of this number for various parameter regime of 
the Bose-Hubbard model. 

To obtain a theoretical understanding of the dynam- 
ics of the Bose-Hubbard model when the interaction pa- 
rameter is quenched from C/j to Uf, we use exact di- 
agonalization to obtain the ground state of the system 
for U — Ui and all states and energies of the system at 
U = Uf. Since the total phonon number of the system 
is conserved, we restrict ourselves to those states in the 
Hilbert space which has a fixed number of total phonons 
(bosons). In the rest of this paper, we take this total 
number of phonons to be equal to the number of sites. 
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FIG. 3: Correlation function of the bosons as a function of 
J/U indicating a finite-sized crossover from the Mott to the 
superfluid phase. 

We denote the initial ground state of the model by |i/>g) 
and the energy eigenstates and eigenvalues at U = Uf 
after the quench by |a) and E a respectively. The wave- 
function of the system at a time t after the quench at 
t = can then be found by solving the Schrodinger 
equation iHdt\ip(t)) — H\tp(t)) with the initial condi- 
tion \ip(t = 0)) = \ipc)- Expanding the wavefunction 

|"0(*)> = La C aWI")> 0ne nnci S 

W)) = £c°e-^|a), C °Hah/> G ). (3) 

a 

The time evolution of any operator Oi, where i is the 
site index, can be obtained from this wavefunction as 
(Oi)(t) = T,a,p c l* c y {Ea ~ Ef>)t,h { a \Oi\P)- In the rest 
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FIG. 4: (Color online) Left panel: Count rate of photons detected per second as a function of time (in units of on-site interaction 
strength). Right panel: Maximum variation of count rate of photons detected per sec. at different initial values of J/U before 
the applied quench. 



of this work, we shall focus on the time evolution of the 
projection of the boson number density on the motional 
ground state of the phonons which corresponds to Oi = 
hoi and is given by 



n i0 (t) 



a,f),k 



C4cos[(E a - E )t/h]dl, ko dl o0 , (4) 



Where d ak = 

ergy eigenstate 



i a \^k ) 1S the projection of the en- 
a) onto the state \^\ ) which has zero 
phonons at the i th site, and the sum over ko indicates a 
sum over all such states in the restricted Hilbert space. 

The equilibrium phase diagram of the Hubbard model 
as a function of J/U, as obtained by exact diagonaliza- 
tion is shown in Fig. [3] where the order parameter cor- 
relation Aij = (b\bj) is plotted as a function of J/U . 
We find that the plot demonstrates a gradual finite-size 
crossover from the Mott phase (A = 0) at small J/U to a 
superfluid phase as J/U is increased. To study the effect 
of local quench, we now prepare the system in its ground 
states with several different U, switch U — > —U at the 
i th site with the peizo actuator, and measure the photon 
count Ri(t) which reflects n^t) as a function of time 
(Eq. [2]). The resulting photon number fluctuation as a 
function of time after the applied quench is shown in the 
left panel of Fig. 2J We find, in accordance with earlier 
theoretical predictions in Rcfs. 



13, list, that the ampli- 



tude of oscillations of the phonon distribution is maximal 
when U is in the crossover region as shown in the right 
panel of Fig. [4] Further, we find that the response of 
the system to the quench is minimal in the Mott phase 
since rii approximately commutes with H for small J/U 
(the commutation is exact for J = 0) and hence can not 
change appreciably as a function of time. From Fig- El we 



estimate that it is feasible to detect a variation of about 
30% in the photon count rate (in a count rate of about 
24, 000 photons per sec) for different interaction times in 
the crossover regime [1_9| . This is well within the present 
experimental reach. 

In conclusion, we have suggested an experimental pro- 
posal of detecting time evolution of the phonons in a 
motional ground state of phonons in an ion trap system 
subsequent to the local change in the effective phonon 
interaction in one of the ion sites. Our work shows that 
such an experiment would shed light on the properties of 
the local dynamics of the Bose- Hubbard model obeyed by 
these phonons which has not been experimentally stud- 
ied before. Our supporting theoretical analysis can be 
easily generalized to initial thermal states of the bosons, 
to quenches with finite ramp rates, and to other models 
emulated by these systems. 
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